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In collisional fluids, a number of key processes rely on the frequency of binary collisions. 
Collisions seem necessary to generate a shock wave when two fluids collide fast enough, to 
fulfill the Rankine-Hugoniot relations, to establish an equation of state or a Maxwellian 
distribution. Yet, these seemingly collisional features are routinely either observed or as¬ 
sumed, in relation with collisionless astrophysical plasmas. This article will review our 
current answers to the following questions: How do colliding collisionless plasmas end-up 
generating a shock as if they were fluids? To which extent are the Rankine-Hugoniot re¬ 
lations fulfilled in this case? Do collisionless shocks propagate like fluid ones? Can we use 
an equation of state to describe collisionless plasmas, like MHD codes for astrophysics 
do? Why are Maxwellian distributions ubiquitous in Particle-In-Cell simulations of col¬ 
lisionless shocks? Time and length scales defining the border between the collisional and 
the collisionless behavior will be given when relevant. In general, when the time and 
length scales involved in the collisionless processes responsible for the fluid-like behavior 
may be neglected, the system may be treated like a fluid. 


1. Introduction 

Astrophysical plasmas are usually collisionless. It means that the mean free path for 
binary Coulomb collisions is usually much larger than the typical size of the system. 
Another way of putting it consists is stating that the average time between electronic 
Coulomb coll isions oc NjyU!~^, where Njj is the number of electrons per Debye sphere 
( Spitzeil 2013 ). is much larger than the other time scales of the problem. In collisionless 
plasmas with No oo, the Spitzer collision frequency goes to zero, and many phenomena 
can be studied by completely neglecting it. 

Another important consequence of collisionless-ness is that collisions can no longer 
thermalize the system on short time scales. Because binary collisi ons are responsible for 
the convergence of the distribution function toward a Maxwellian ( Chandrasekhar] [l 943ll 
on a “Spitzer time scale”, the absence of collisions suppresses the most common ther- 
malization agent. It may therefore be problematic to use an equation of state for the 
medium, since the very existence of such an equation assumes at least a local thermal 
equilibrium. Yet, astrophysicists make an intensive use of MHD codes requirin g an equa¬ 
tion of state ( Hawley et al. 1984 ; Stone fc Normanlll992 ; Sadowski et al. 2014 ). To which 
extent is such a strategy appropriate? Although we don’t have any definite answer so 
far, progresses have been made in recent years pointing toward faster-than-Spitzer ther- 
malization mechanisms, which could justify treating collisionless plasmas like collisional 
fluids in appropriate settings. 
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Collisionless shocks illustrate also the fluid-like behavior of collisionless plasmas. When 
two collisional fluids collide fast enough, or when a large amplitude sound wave is 
launched in one fluid, a shock wave is generated ( ZePdovich fc Raizei 20021) . At the 
fluid, macroscopic, level, this wave represents a traveling discontinuity. And macroscopic 
quantities like density, temperature or pressure on both sides of the discontinuity obey 
the so-called Rankine-Hugoniot (RH) relations. It turns out that the encounter of two 
collisionless plasmas also results in the formation of a shock, although both plasmas never 
“hit” each other, but start passing through each other instead. Furthermore, simulations 
show a fast th ermalization (or a t least a thermal component ) of the medium when cross - 
ing the front ( Soitkovskv 2008 : iLaoenta et al. 20071120091 : ISironi fc Soitkovsl^ 2011 ). 
The RH conditions a re satisfied to a rea sonable accuracy, whether for simulations or for 
the earth bow shock ( Farris et al. 19921) . 


We thus find several instances of fluid-like behavior of collisionless plasmas. Whether 
it be the issue of thermalization or the shock problem, it seems that some collisionless 
mechanisms can achieve what binary collisions cannot. We will here review the current 
status of our knowledge in these matters. 


2. Fluid-like shock physics 

When a shock propagates, upstream particles need to dissipate part of their kinetic 
energy for the fluid to slow down in the downstream. For a collisional fluid, this is achieved 
through an increase of the numb er of inter-particle collision s. This is why the shock front 
is a few mean-free-paths thick (jZerdovich fc Raizen l2002i) . that is, a discontinuity at 
the fluid level. Let us now turn to the earth bow sh ock. Here, we know from in situ 
measurements that its front is about 100 km thick ( Bale et al. 20031 : ISchwartz et al. 


20111) . But the proton mean-free-p ath at the same location is comparable to the sun- 


earth distance (iKasoer et aLN2008l) . Or turn to the supernova remnant SN 1006. Here, 
the s hock front is about 0.04 pc (1.2 x 10^^ cm ) while the mean-free-path is 13 pc (4 x 10^® 


cm) ( Bamba et al. 2003 : Takabe et al. 2008h . 

It is thus clear that the mechanism ensuring dissipation at the shock front cannot come 
from binary collisions for these astrophysical shocks. I ndeed, there was a time when the 
mere existence of such shocks was simply questioned (Sagdeey 1966t Sagdeev fc Kennel 
199lh . We shall now see how they are formed, and how they result in structures fulfilling 
the RH jump conditions. 

Although this rev i ew is n ot about collisionless s hocks per see (see iTreumannI ( 20091) : 
Bykov fc TreumannI ( 201ll) : iBalogh fc Treumann ( 2013 ) for such reviews), we briefly 
mention in Appendix the various kind of collisionless shocks one can encounter. 

It is worth specifying that the fluid-like behaviors we’re commenting relate only to 
unmagnetized shocks in pair and electron/ion plasmas. Yet, to our knowledge, no simu¬ 
lation of colliding plasmas that would produc e a shock for collisional fluids, has failed so 
far to produce one with collisionless plasmas ( Sironil 2014 ). 

We now turn to our first example of fluid-like behavior. How does the encounter of two 
unmagnetized, collisionless plasmas, result in the formation of a shock, like it happens 
for fluids? The detailed theory of collisionless shock formation we are about to review for 
two simple cases will help understand. 


2.1. Unmagnetized pair plasmas 

The setup consisting of two identical colliding pair plasmas is attractive for several rea¬ 


sons: 
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(a) 

(b) 


(c) 


(d) 


Shock Front Shock Front 

Figure 1. Formation of a collisionless shock via the encounter of 2 pair plasmas, (a) The 2 plas¬ 
mas shells are heading toward each other, (b) They start interpenetrating, and the overlapping 
region turns unstable, (c) Instability saturates, generating a strong enough turbulence to block 
the incoming flow. The density in the central region increases until RH conditions are fulfilled, 
(d) Two shocks form and propagate. For simplicity, identical plasmas are shown. 



Turbulence, density increases 


(a) The absence of particles of different mass renders both theory and simulation much 
simpler. 

(&) B y virtue of this mass issue, pair plasmas don’t exhibit a Debye sheath at their 
border ( Gurnett fc Bhattacharie3 200^. In case two of them collide, they will always 
have enough energy to start interpenetrating. 


(c) Pair plasmas are directly relevant to some astrophysical settings (lUsovI 11992 : 
BrodCTic^_et_a||2012 ^ Sironi fc Giannio^l2014 ), and available in the laboratory ( Chen et al. 
2 OI 0 I : Sarri et alT 2013h . In addition, experiments mimicking p air plasmas can be per¬ 
formed with positively a nd negatively charged Cqq molecules (lOohara fc Hatakevama 
20031: lOohara fit r7,/.ll2nn,^ . 

In the absence of binary collisions, two identical plasmas heading toward each other at 
velocity ±vo don’t collide (see hgure[T]). Instead, they simply pass through each other. 
How is it then that a shock eventually forms? The overlapping region could grow indef¬ 
initely, atjeas^_imdl_the_^itzer time, if it were not for counter-streaming instabilities 
I Kahn 1958; Bunemaiil 1964h . Being formed by two counter-streaming plasmas, instabili¬ 
ties will grow in this region, saturate, and generate a turbulence. Though many instabil¬ 
ities do grow, as explained in Appendix lA.il the Weibel instability is the fastest growing 
one in these circumstances, which is why these shocks are called “Weibel mediated”. 

If the randomization length of the incoming flow L^, is smaller than the size of the 
central regio n at saturation Ls, the plasma which keeps entering the overlapping region 
stops there (|Bret et a/.l 1201311201411 . As a result, the bulk velocity in the central region 


IS now zero. 

At this stage, the system presents therefore a macroscopic discontinuity in velocity 
space, but the central density is still the sum of the densities of the two shells. The time 
needed to get to this stage from the first contact of the plasmas is given by. 


= Cn, 


( 2 . 1 ) 


where Sm is the growth rate of the dominant instability (see Appendix lA.ll) and H the 
number of exponentiations until saturation. Provided the soon-to-form shock will have 
a density jump T given by the RH relations (see section |3]), the density jump will now 
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increase from 2 to F. If it took a time Tb to bring enough material in the central region 
for a density jump 2 to build-up, the time needed to reach F, namely the shock formation 
time, will be 

Tj ^T, + (F-2)r, = (F-1 )t,. (2.2) 

Therefore, provided one is interested in time scales longer than Tf and length scales 
larger than Lg ^ ^voTg, we recover a fluid-like behavior: the two shells “collide”, and two 
opposite shocks are generated in each one. 

2.2. Unmagnetized electron/ion plasmas 

Things are different when particles of various masses are involved, especially by virtue 
of the large proton-to-electron mass ratio. On the one hand, various instabilities may 
develop on their proper time scale, so that the building up of a turbulent region is more 
involved than for pair plasmas. On the other hand, the Debye sheaths bordering these 
plasmas can interfere at low energies, and alter the shock formation. 

For low impacting energies, i.e., if the initial kinetic energy of our identical shells is not 
high enough for protons to overcome the pote ntial barrier of the other shell, a collisionless 
electrostatic shock is form ed within a time ( Forshmd fc ShoiiS 1971)1: Dieckmann et a,I 
l201,ll: IStockem et al. 2014 1. 


Tf ^ lOyo^^w-i. (2.3) 

Note that the impacting velocity may be relativistic for high-enough temperatures, hence 
the Lorentz factor 70 in the equation above. Also, the dominant time-scale is now the 
protonic one, which is why the formation time scales like the proton plasma frequency 
ujpi- Protons also govern the width of the shock front, which is eventually a few c/ujpi 
thick. 

As the impacting energy increases, a transition region h as to be crossed whic h is not 
straightforward to describe. Here, we can refer the reader to lStockem et al\ (1201411 . Then, 
for high enough energies, roughly (70 — l)meC^ ^ kgTg, we recover a situation similar to 
the pair plasmas case. Because there is enough energy to overcome the potential barriers, 
the two plasma shells start to overlap. The counter-streaming electrons first turn Weibel 
unstable, and their dominant instability saturates. The counter-streaming protons then 
turn Weibel unstable, and the relevant sy stem at this stage consists in two opposit e 
proton streams over a bath of hot electrons ( Shaisultanov et al.ll2012l : IPavis et aLll2013l l. 

Denoting again dm the growth rate of the dominant instability of this system, we can 
repeat the reasoning done previously for pair plasmas, and derive an equation similar to 
( 1 ^ for the formation time. Note however that some theoretical work is still needed to 
fully explain the simulations. One reason is that the growth rate 6m depends on the tem¬ 
perature of the hot electrons bath. It appears from the simulatio ns that they are heated to 


a temperature much higher than their i nitial streaming energy (IGedalin. M. et al.N2012 


Plotnikov et al. 2013 : Davis et'^l2013ll . Although various theories have been proposed 


to explain this anomalous electron heating, the situation is not clear yet. 

To conclude this section on shock formation, we can say that a collisionless shock 
forms when two collisionless plasmas collide, over time and length scales given by the 
equations above. We therefore recover a fluid-like behavior, though the underlying physics 
completely differs from the collisional case. 


2.3. Shock stability and stationarity 

Once the shock has formed, the RH conditions determine its velocity, both in the fluid 
and the collisionless case (see section [3|). A key issue which has been established long 
ago in the fluid case is the one of the stability of the shock front. Whether we deal with 
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Figure 2. Non-stationarity of a collisionless shock. The front propagates from left to right. A 
is the numerical grid size. The evolution of the field at the front evidences the non-stationarity 
of the shock. The upper -left green rectangle pertains to the moving piston generating the shock. 
From l^ng et a[\ ll2009li . 


the non-relativistic ( D’yakov 1954 Kontoroyich 1958; Wouchuk fc Lopez Cavadal 2004ll 
or the relativistic case ( Anile fc Russolll987l) . a corrugation applied to the shock front 
in an ideal fluid does not grow. Though instabilities have been evidenced for shocks i n 
substances described by a van der Waals equation of state ( Bates fc Montgomery 2000ll . 
shocks in ideal fluids are stable. 

What about the collisionless case? To our knowledge, no theoretical work has proved so 
far the front stability or instability against all kinds of perturbations. The challenge here 
is that electromagnetic perturbations have to be considered, in addition to the shock 
front corruga tion. Regarding stabil ity against electromagnetic fluctuations for electro¬ 
static shocks. [Stockem et al. (2014) could relate it to the electrostatic/electromagnetic 
shock transition. 

The process of collisionless shock re-formation could indicate a departure from the 
fluid stability, or at least station arity. This phenomenon was f irst noticed in ID simula- 
tions ( Lembege fc Dawsoi]ll987 b before it was observed in 2D ( Lembege fc Savoi'nilll992 ) 
and in nature on t he earth bow shock ( Mazelle et 20inl l. Also referred to as shock 
“non-stationarity” ( Lembege et aLlliooi b this process occurs for collisionless shocks prop¬ 
agating quasi-perpendicularly to an ambient magnetic field. 

This self-reformation is due to reflected ions which accumulate in front of the shock. 
It results in a non-stationary evolution of the shock front on a time scale given by the 
upstream ions gyrofrequency ^ and a length scale about the ions gyroradius A^, as 
evidenced on figure [2l For this figure, the shock has been generated by a piston moving 
rightward, and visible on the plot through the upper-left green rectangle. Like fluid shocks 
( Landau fc LifshiQl201,8 b collisionless shocks can also be generated by a moving piston. 

Non-stationarity is not the fruit of some external factor. It is part of the proper shock 
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dynamics. Provided one is interested in long enough time scales and large 

enough length scales (s> Xi), a. non-stationary shock behaves like a fluid one. So far, it 
has been found stationary. 


2.4. Downstream ion-electron thermal equilibration 

The issue of ion-electron thermal equilibration in the downstream of some collisionless 
shocks could equally imply some departure from the collisional behavioiQ. Based on 
observations of some Supernova remnan t s, gal axy cluster shocks and the terrestrial and 
Saturnian bow shocks. iGhavamian et al. ( 2013h found that the ratio T^/Ti of the electron 
to the ion temperatures falls like A4~^, where A4 is the Mach number. An analytical two 
fluids model established Tg/T^ ~ 1 f or Af 5 . 2 , T^ /Ti melrrn for A4 > 60, and T^/Ti oc 


A4 ^ in the intermediate regime (|Vink et al\ 120141) . Little is known howe ver on the 


micro physics involved in equilibrating the ions and electrons temperatures (|Matsukivo 
I 2 OIOII . Time will therefore be needed to draw definite conclusions on the collisional status 
of this feature. At any rate, the ultimate collisionless picture could be equivalent to a 
two-fluids model, instead of a single-fluid one. 


3. Rankine-Hugoniot relations and eqnation of state 

The foundations of the RH jump conditions reside in the conservation of matter, mo¬ 
mentum and energy. For the fluid case, they are si mply derived writing that t hese quan¬ 
tities are conserved when crossing the shock front (jZel’dovich fc Raizerll2002l) . To which 
extent can they remain valid for a collisionless shock? 

Since matter carries all conserved quantities (we exclude radiative shocks), we need 
to make sure all the matter upstream ends up downstream in the collisionless case. This 
is verified, but not perfectly because some particles from the upstream can be reflected 
by the shock front, and come back upstream. Also, other particles may be accelerated 
near the front, escaping the system nearly isotropically. A numerical i nvestigation of 


these effects demonstrated a departure from RH of a few percent (see IStockem et al 
( 2012h and concluding remarks). For shocks which are efficient particle accelerator^. 


up to 10 %- 20 % of the upstre am kinetic energy can be converted to energetic particles 
( Caprioli fc Spitkovskvl 120141 . As a result, the downstream is found 20% cooler than it 
would be for a fluid (see figure [3]). 

We can therefore state that up to a ~ 10 — 20% accuracy, all the matter upstream 
of a collisionless shock passes downstream, together with the momentum and energy it 
carries. Yet, the RH conditions under their usual form need more than the conservation 
laws. They need an equation of state of the form, 


p = (7 - m 


(3.1) 


where P is the pressure, 7 is the adiabatic index and U the energy density of the medium. 
As showed in Appendix m such an equation of state does not require a Maxwellian 
distribution. An isotropic distribution function, in the reference frame where the plasma is 
a rest, is enougl0. In simulations, the upstream distribution is an input usually considered 
isotropic. Regarding the downstream distribution, turbulent magn etic fields are presen t 
both in the downstream of electromagnetic and electrostatic shocks ( Stockem et al.l[2014 ). 


t See also l|R.avmond et aZ.ir2010l '). 

j In the magnetized regime, it seems these are the quasi-parallel shocks, whe re the angle 6b 
between the shock f ront normal and the external field is not “to o close” to 'k 12 llLemoine et aZ.I 
I 2 OO 6 I : [Capriolil[2014l l. For non-relativistic shocks. FCapriolil (l2014l 'l gives 9b iS 7r/4. 

^ I thank Guy Pelletier for pointing this out to me. 






































Figure 3. Downstream ion energy distribution in terms of time in a shock simulation. A 
Maxwellian is clearly evidenced, together with a power-law tail extending to larger energies 
at later times. The spectrum is multiplied by to emphasize the agreement with the scaling 
predicted by the acceleration mechanism. Time is measured in unit of the ion cyclotron 
frequency in the external magnetic field. The energy E is normalized to E sh ~ , where Vbh 

is the upstream fluid velocity in the downstream reference frame. From [Caprioli fc Spitkovskvl 
(1201411 . 


that are capable of nearly isotropizing the distributions. It is therefore not surprising to 
observe the near fulfillment of the RH conditions in simulatio ns of collisionless shocks 
(see for example Soitkovs^ ( 2008 1: Sironi fc Soitkovskv ( 2009ll l. 

Yet, it seems we not only observe isotropic distributions. We observe Maxwellian dis¬ 
tributions. Some processes must therefore be at work to operate such a collisionless 
thermalization. This is the topic of the next section. 


4. Collisionless thermalization 

Figure |3] displays the time evolution of the downstream ion energy distribution in a 
shock simulation. A remarkably stable Maxwellian component is clearly visible at low 
energies. At t = 2500a;“^, its temperature is 80% of what it would be for a fluid shock. 
Anoth er component, a power-law with a cut-off, is present as a result of p article acceler¬ 
ation ( Fermi 19491 : Krymskii 1977 : Blandford fc Ostrikeilll978 : Bell 19781) . 

Where does this Maxwellian come from? As already mentioned, the magnetic fields 
which are present either in the shock front or in the near downstream, can explain 
isotropization, but thermalization, i.e. Maxwellianization, is achieved over time scales far 
shorter than the Spitzer one. In the absence of collisions, the distribution function is 
solution of the Vlasov equation, 


dt 


+ V 


dr 


q 


et-vxbV 1^ = 0, 

C / 


(4.1) 


which does not demand a Maxwellian. What is then its origin? 

It is remarkable to notice that the same enigma equally arises when studying the dis¬ 
tribution of stars’ velocity in a galaxy. As noted bv lLvnden-Bell ( 1967t ). “The remarkable 
regularity in the light distribution in elliptical galaxies suggests that they have reached 
some form of natural equilibrium. However, estimates of the normal star-star relaxation 
show that it is too weak to establish equilibrium in the time available”. The concept 
of “violent rela xation” introduced by Lynden-Bell to explain this phenomena, is still 
debated today ( Arad fc Lvnden-Belll 20014 Binnev fc Tremaine 201lh . 
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Coming back to plasma physics, the anomalous occurrence of Maxwellian distribution 
functions has been known for nearly one century as “Langmuir’s Paradox”. In 1925, 
Irving Langmuir observed the persistence of a Maxwellian distribution in an experiment 
where “the number of collisions of the electrons with each other and with atoms were 
far too few to maintain the observed distributi on” (iLangmuid 192511 . The solution to this 
problem would require decades of efforts, and Gabor et all ( 1955 1 even referred to it as 
“one of the worst discrepancies known in science”0- 
Kadomtsev (196^ discussed this process in terms of turbulent plasma heating. This 


author invoked the strong turbulence generated by beam plasma instabilities to explain 
collisionless thermalization. Yet, no definite theory could be laid out at that time, as 
“the real picture of the excitation of waves by powerful beam and the interaction of the 
waves with the particles and with one another may turn out to be much more complete, 
and to require a c onsiderable development of the theory for its complete description” 
( Kadomtsev! ( 1965ll . p. 131). 

Without developing a detailed theory. iHovauxI ( 1968ll noted that the randomizing ef¬ 
fect of plasma turbulence could be modeled considering particles are subject to a macro¬ 
scopic force and a microscopic, stochastic, one. The resulting Langevin equation models 
a Brownian motion in velocity space which provi des a Maxwellian distribution, provided 
the s tochastic force has a Gaussian distribution ( Chandrasekhad 19431 : Dunkel fc Hanggi 
200^. Such is the way Dieckmann et al\ ( 2006ll give account of the Maxwellians found 


in the beam-plasma simulation they study. 

After more than 80 ye ars of research, a breakthrough has recently been made b;v 


Alter more tnan oU ye ars ol researcn, a Dreaxtnrougn nas recentl y been made by 
Baalrud_et_^ ( 20091 2010ll . Starting either from the BBGKY hierarchy ICercignanilll97CT 


IchimanJ 1973ll or from the Klimentovich equation ( Klimontovich 1967f) . these authors 


derived an “instability-enhanced” collision term for Eq. which accounts for wave- 
particle scattering due to plasma instabilities. This non-relativistic theory neglects the 
magnetic field produced by the charges in motion. When applied to the conditions of 
Langmuir’s experiment, the enhancement of electron-electron collisions produced by the 
ion-acoustic instabilities produces a convergence toward “nearly a Maxwellian” (more on 
this in the conclusion) 100 times faster than the Spitzer time, explaining the paradox. 
The number 100 depends on various factors, like the plasma density and temperature, 
the experimental setup and the kind of unstable modes excited. 

Noteworthily, convergence to a Maxwellian requires a collision operator with both 
a drag-term an d a diffusion-term. The quasi-linear theory like the one developed by 
IVedenovI ( 1963h only displays a diffusion term. It misses a drag term to describe ade¬ 
quately the convergence to a stationary distribution. Regarding the Fokker-Planck equa¬ 
tion, it is known it reduces to the quasi-linear theory when the two t erms of its collision 
operator are computed neglect ing the recoil of the particles’ scatterers (|Blandford fc Eichlei 
1987 : Peeters fc Strintzill2008l l. As such, it cannot describe either the convergence toward 


a stationary distribution. 

At any rate, an instability generated turbulence seems required to achieve fast ther¬ 
malization. This is related to the fact that the instability-enhanced collision operator 
fulfills the "H-theorem, and vanishes for a definite distribution. Let us remind that in 
contrast, a stable Vlasov plasma, w hen perturbed even non-linearly, evolves at constant 
entropy (IPitaevskii fc Lifshitzlll981ll and comes back to a different distribution function 


f See iBaalrudI (|2ni(lf) for a detailed history of the Langmuir’s paradox. 
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than the one it had in the first place ( Villanil[2014l l. Simply put, Landau damping does 
not restore a MaxwelliailJ. 

Turbulent settings, like the ones arising from the development of plasma instabilities, 
favor the collisional behavior of collisionless plasmas. One the one hand, turbulence leads 
to an isotropic medium. On the other hand, plasma turbulence is l ikely to produce the 
kind of enhanced collisionality described by 
non-relativistic, electrostatic case. 


Baalrud et al. ( 20091 2010l l. so far in the 


5. Instabilities and related phenomena 

We now comment on parallels between some instabilities and related phenomena in the 
collisional and collisionless regime. Clearly, in order to compare any collisionless behavior 
to some fluid counterpart, the instability under scrutiny must exist in the fluid world. 
It means it must exist in a system which can be modelled by a single fluid at a given 
position in space. Excluded then are the plasma streaming instabilities, like two-stream 
or filamentation, which by definition involve the counter-motion of several fluids at the 
same place. Note than these can be described quite successfully by a two-fluids model 
(jBret fc Deutschl 120061) . Yet, we here restrict to the parallel with a single fluid. 

In the case of the fluid Kelvin-Helmholtz instability, a corrugation at the interface be¬ 
tween two shear flows grows as a result of vorticity dynamics on each side of the distur¬ 
bance ( Batchel^l200(]ll . For flows of the same density with velocities ± V, the fluid growth 
rate for disturbances with wave vector k || ±V is Skf.f = ^kV ( Landau fc Lifshit^ 
2 OI 3 II . The same instability has been studied both theoretically and numerically for 


two semi-infinite colli s ionless plasmas d rifting in opposite directions ( Griizino'^ I 2 OO 8 I : 
Grismaver et all l201,ll : lAlves l liHl . For k <C ujpfV, the growth rate is simply 


‘2dKF,F = kV. But here, the underlying mechanis m is similar to the tw o-stream instabil¬ 
ity, which has also S = kV for large wavelengths ( Mikhailovskiil [l973l . The coincidence 
of the growth rates is striking, in spite of apparently different instability mechanisms. It 
is probable that some deeper physical connection can be found between the fluid and the 
collisionless plasma KH instabilities. 

Instabilities have been also considered in connection to the collisional behavior of 
ions in incompressible, collisionless turbulence. Again within shear flows, here magne- 
tized, a turbulence is generated out of the growth of the firehose and mirror instabilities 
( Kunz et al. 20 l3l. or out of the growth of the mirror and the ion cyclotron instabili¬ 
ties (|Riauelme et A collisional behavior is retrieved. For example, whether one 

deals with the firehose or the mirror instability, the instability itself generates an effective 
scattering rate ma intaining the pressure anisotropy to the level of a collisional plasma 
( Kunz et al. I [2^ . 

The magnetorotational instability is of great relevance for the physics of accretio n 
disks. It has been studied first within the MHD framework ( Balbus fc Hnwle^ 1998b 


PIC simulations evidenced a behavior very similar to MHD simulations (IRiauelme et al 


l2012h . Theoretical works found identical kinetic and MHD stab ility criteria, while th e 
growth-rates agreement depends on the plasma magnetization (IQuataert et aZ.M2002ll . 
Indeed, the MHD-like behavior of the MRI is tightly related to the presence of kinetic 
plasma instabilities th at suppress the gro wth of pressure anisotropies, thus mimicking 
the effect of collisions ( Sharma et al. 2006[l . 

Also relevant to the physics of accretion disks is the problem of ion-electron tem- 


1: Which is the reason behind plasma echo (iGould et al\ (Il967f) . iPitaevskii fc Lifshitd (Il98lfl . 
p.141). 
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perature e puilibration. In low-lum inosity accretion disks, the ions are hotter than the 
electrons ( Yuan fc Nar avail 2014l i. In such collisionless environments, binary collisions 
could still equilibrate th e temperatures , but on the extra-long Spitzer time scale. Hence, 


Sironi fc Naravan (2014); Sironi ( 2014h investigated how ion driven plasmas instabilities 


can transfer energy to the electrons. In this respect, it seems there is “solid evidence” 
that the ion cyclotron instability tend to equilibrate ionic and electronic temperatures, 
if the latter are initially much cooler than the former (Tpe < 0.2Toi). Further works will 
be needed to assess to full parameter range, allowing for example to introduce an extra 
electron heating term in MHD codes. 


6 . Conclusion 

Collisionless plasmas behave quite like collisional fluids in some conditions. Like in 
fluids, shocks can be generated launching them against each other. To the extent that 
the downstream distribution function is reasonably isotropic, an equation of state can be 
used, and RH jump conditions applied with good accuracy. As it propagates, the shock 
may reform and/or accelerate particles, so that its stationarity is only achieved up to 
certain time and length scales. 

The progressive building-u p of the high energ y power-law tail observed on Fig. [3] 
eventually reaches satu ration ( Caprioli et al\ \20l4) . For weakly magnetized shocks with 
a < 10“^ (see Eo. lA ip . ISironi et al. ( 20131) found that the maximum energy of the accel¬ 
erated particles grows like before it saturates when these particles’ energy reaches 
about 10% of the flow. The 20% departure from the downstream RH temperature could 
thus be an upper-bound. Indeed, some Supernova Remnants have been aro und for mil- 
lenni|0, su ggesting they could h ave reached a kind of steady state. Even if ISironi et al 


([20131) and lCaorioli et ali (|2014r) show how this can happen for the shock they study, the 
same conclusion is yet to reached for all kinds of collisionless shocks. 

Also, the “faster-than-Spitzer” Maxwellianization observed in experiments or simu¬ 
lations has been explained, so far for the non-relativistic electrostatic case only, as an 
instability-enhanced collisional effect. The system is still collisionless, but particles scat¬ 
tering on an instability generated turbulence, drive a fast thermalization. 

In spite of these numerous fluid-like features, some words of caution seem appropriate: 

• Using an equation of state requires an isotropic distribution (see Appendix IbI) . 
Yet, such isotropy can be jeopardized in strongly magnetized settings. In this respect. 
Chew et~ai] ( 19561 ) tried to derive the MHD equations, starting from the Vlasov equation 
with a large Lorentz force. Their attempt was dimmed “not entirely successful” due to 
the anisotropy brought by the magnetic field lines. 

• Even for the instability-en hanced thermalizatio n, the distribution function converges 
only “nearly to a Maxwellian” ( Baalrud et al. 20101) . because the collision operator van¬ 
ishes for distributions only close to a Maxwellian. This could be the reason why an 
electrostatic shock produces a flat-top electron distribution in ID. In ID-electrostatic 
the distribution is constrained by the fact that phase space paths cannot intersect. A 
2D simulation with the same initial conditions generates ion acoustic turbulence down¬ 
stream of the shock and that one prod uces rapidly an electron distribution with a sin gle 
maximum that looks Maxwellian-like ( Dieckmann et ai]l2013l:lDieckmann et a/. l2013l ). 

Instability-enhanced collision frequency seems thus the key to fast thermalization. Now 
that Langmuir’s paradox has been solyed, one can wish further works will characterize the 


t SN 1006, already mentioned, was born in AD 1006 (iM urdin fc MurdinI 19851) . The Supernova 
Remnant GSH 1380194 is 4.3 million years old (|Stil fc Irwinll200ll ). 
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Figure 4. Zoology of collisionless shocks. Different kind of shocks develop according to the 
magnetization, and the height of the Debye sheaths bordering the colliding plasmas. This clas¬ 
sification could be enriched considering even more ingredients (see text). 


turbulence needed to generate Maxwellian distributions, beyond the initial Langmuir’s 
experimenta l setup. Since (quasi) therm alization is observed in PIC simulations beyond 
the scope of Baalrud et all (|200i] , 20ir)[l , it could be possible to extend the proof to the 
relativistic regime, or to the case of a full electromagnetic turbulence. 

Beyond this, the similarities between the fluid and the collisionless Kelvin-Helmholtz 
instabilities allow to think that some general physical principles could be at work behind 
these apparently different processes. On the long term, maybe such principles will allow 
to draw a clear quantitative line between the collisionless and the fluid behaviors of 
diluted plasmas. 


7. Acknowledgements 

The author acknowledges the support from grant ENE2013-45661-C2-1-P from the 
Ministerio de Educacion y Ciencia, Spain. Thanks are due to Didier Benisti, Eric Black¬ 
man, Damiano Caprioli, Mark Dieckmann, Laurent Gremillet, Thomas Grismayer, Bertrand 
Lembege, Ramesh Narayan, Guy Pelletier, Aleksander Sadowski and Gustavo Wouchuk 
for inputs and enriching discussions. 


Appendix A. Shock zoology 

Because plasmas are composed of charged particles, they can interact with electro¬ 
magnetic fields. This feature, combined with the Debye sheath likely to border a plasma, 
introduces a rich variety of cases when considering the encounter of two plasmas. The 
typology related to shocks generated by colliding 2 plasmas is schematically represented 
on figure m _ 


• T he first important parameter is the magnetization. At low magnetization (jNiemiec et al. 
2 m 2 I) . namely 

^ — / 1 \ 2 ^ 

(70 - IjnoTOC^ 

randomization of the upstream at the front is provided by the magnetic field generated 






















12 


A. Bret 


by the Weibel instability (see section EH. For high magnetization, it rather comes 
from reflection on the compressed field at the front. In this respect, besides the field 
comp ression due to the RH conditions including the magnetic energy (jTidman fc Krall 
197lh . so me plasma instabilities can provide an additional enhancement of the field at 


the front ( Bell 2004[) . 

• Then comes the height of the Debye sheath, compared to the kinetic energy. Because 
ions are much heavier than electrons, the electronic density does not drop to zero like the 
ionic one at the border of an electron/ion plasmas. Such plasmas are instead bordered by 
a Debye sheath displaying a potential ju mp kj^Tp/q high, where Tp in the electronic 
temperature and q the elementary charge ( Gurnett fc Bhattacharie3l2005ll . If the initial 
kinetic energy of the colliding plasmas is high enough, that is, (70 — l)nomc^ ksTe, 
both interpenetrate and an electromagnetic shock will be generated. On the contrary, 
if the interaction is mediated by the one of the Debye sheaths, an electrostatic shock 
is formed. Although not formally proven, this is probably the reason why electrostatic 
shocks cannot exist in pair plasmas with equ al temperatures for electrons and positrons 
( Oohara fc Hatakevamal[2003 IVerheest 2005ll . But if electrons and positrons have differ¬ 
ent temperatures, a Debye sheath can fo rm, which is probably connected to the possibility 
of electrostatic shocks in such plasmas (iDubinov et aLll2006llii . 

Although their formation mechanism slightly differ (see sections 12.11 & 12 . 21 ) , shocks in 
pair or electron/ion plasmas can therefore be classified according to the same typology 
when focusing on the Debye sheath potential. Still, figured falls short of exhausting every 
possibility. To do so, we would have t o evoke the encoun t er of t wo different plasmas (pair 
vs. electron/ion, different densities... Dieckmann et ^ ( 2014l l'l. We could also contem¬ 
plate plasma shells inclu ding more than 2 species, like electron -|- positron & proton, or 
electron -|- ion & proton ( Dieckmann fc Bi^ 2009ll . 

For magnetized shocks, yet another key parameter is the angle 9b between the shock 
front normal and the external field. For example, non-relativistic perpendicular shocks 
{9b = 7 r/ 2 ) exhibit interesting non-stationary features that are discussed in section [2.31 
Also ultra-relativisti c perpendicular shock s are not good Fermi accelerators, while parallel 
shocks {9b =0) are ( Lemoine et al. 200611 . Finally, cases are when the energy of radiation 
is large enough to significantly enter the energy budget and alter the jump conditions. 
Such “radiative” shocks, collisionless or not, a re encountered in astrophysical conditions 


and have equally been studied in literature (|Zerdovich fc Raized l2002t iBouauet et al 

2n0(t . 


A.l. Dominant instability and turbulence 

We here comment on the dominant instability at work in the shock formation, and on 
the resulting turbulence. Note that we are only interested in these processes to the extent 
that they participate in building up a shock. We thus won’t give much detail on them 
and refer the reader to external references. 

Counter-streaming instabilities are th e reason why the overlapping region does not keep 
expanding, but turns turbulent instead ( Bunem^ll964ll . Counter-streaming systems can 
be destabilized by many kind of instabilit i es, an d the dominant one depends on the 
parameters of the problem (see iBret et 'al,. ( 201fll l for a review). For counter-streaming 
pair plasmas of nearly equal densities, the fastest growing modes are likely to be found 


f Strictly speaking, these references deal with solitary waves, not shock waves. Yet, as nicely 
explained by Sagdeev, both strnctures are closely related since a soliton is a large amplitude 
wave leaving the medium unchanged after it s passage, whi l e a shock d oes affect the medium. 
Simply put, “soliton -|- dissipation = shock” (ISagdee-dllfififil : iKra 
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2008h. These modes are now referred 


with a wave vector normal to the flow (IBret et all _ _ __ 

to as filamentation, or Weibel, modes ( Weibell IQhdT FriedI 1959h . 

In the case of electron/ion plasmas, where the relevant counter-streaming system con¬ 
sists in two opposite protons beams over a b ath of hot electron, Weibe l modes have also 
been found to govern the unstable spectrum (|Shaisultanov et aLll2012h . “Weibel shocks” 
can therefore be found for a wide rang e of plasma types . Yet, it is well kno wn that this in¬ 


stabi l ity is vulnerable to temperature (ISilva et al.ll2002f) or magnetization (jGodfrev et al 


1975t Bret et al. 2006ll . Simulations have been conducted wh ere a flow-aligned magneti c 


field had been setup precisely to stabilize the Weibel modes ( Sironi fc: Soitkovskv 2009ll . 
In this case, other unstable modes than Weibel take the lead. A turbulence is still gen¬ 
erated which results in the formation of a shock. 

Once the shock has been formed, the turbulence behind the front is maintained by 
the incoming upstream flow. This is how the collisionless shock eventually results in a 
self-sustained structure. Pla sma instabilities keep arising near the shock front as a result 


of counter-streaming flows ( Karimabadi et al. 1991 : Kato fc Takabel [20101 : iNakar et al. 


2011 : Lemoine et aLll201^ . But far dow nstream, where the upstream flow cannot arrive. 


the turbulence decays by phase mixing ( Chang et al. 2008t Lemoind 2014 1. 


Appendix B. Basic requirements for an EOS 

We here prove that it only takes an isotropic distribution function to have an equation 
of state like (EH). To this extent, we start performing a ID anal ysis of the pressure in 
a gas with only 1 species. The calculation follows the derivation of iFevnman fc Leighton 
I I977II (see §39-2) who start from the definition of pressure. 

Suppose a ID gas along an x axis, and an obstacle on the axis at a; = 0. Particles 
from its left bounce back against it, pushing it. The gas has the distribution function 
f(x,Vx) (time dependence can be included). When bouncing, a particle with velocity Vx 
transmits the momentum dp = 2mVx to the obstacle. How many of them do so in a time 
dt? All the particles at distance dl from the wall such that dl < Vxdt. The total amount 
of momentum transferred to the wall is thus, 

p<x) pO 

Ap= dvx / dx f{x,Vx) 2mvx, (B1) 

Jo j —Va^dt 

where velocities are integrated from 0 to 00 only, because leftward particles do not hit 
the wall. We assume here that particles hitting the wall are by no means influenced by 
the others, which comes down to the perfect gas hypothesis. 

For dt —>• 0, —Vxdt —>• 0“, and 

/ dx f{x,Vx) 2mvx ^ Vxdt /(O ,Vx) 2mVx, (B 2) 

J — Vxdt 

so that, 


/•oo _ Ap _ 

Ap= dvx Vxdt /(O ,Vx) 2mvx ^ -j- = 2 / dvx Vxf{0 ,Vx) mvx- 
Jo dt Jq 


(B3) 


If f{vx) is isotropic, it is also an even function. Therefore, 2 °° dvx = dvx and 


Ap 

dt 


/ OO 

dvx f{0~,Vx) Vxmvx ■ 

-OO 


Momentum flux at x—0 


(B4) 
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As is well known, the pre ssure is thus related mic roscopically to the momentum flux 
at the left of the obstacle ( Landau fc Lifshi't3l201,'lh . The 3D result is straightforwardly 
derived, and reads (the wall is now T to cc), 


P = y (fvf{0 ,v) (B5) 

If /(v) is isotropic, then 

J (fvfvl = j (fvfVy = J d^vfvl ^ J (fvfvl = ^ J d^vfv'^. (B6) 

Hence, 

P=^ J d^v /(0“,v) =‘^ J d^v /(0”,v) = |t7, (B7) 

'-V-' 

Energy density U 

where the perfect gas hypothesis again is implied, if the energy density is to be solely 
made of kinetic energy. An equation of state of the form 

P = |c/= (7 - 1)C/, (B8) 


is therefore recovered, where 7 is the polytropic index. It is clear that the “2” at the 
numerator comes from the necessity to cancel the ^ factor in the kinetic energy formula, 
and that the “3” at the denominator is related to the number of dimensions of the system. 
Therefore, 7 is still defined through 

7=l + y (B9) 


where d is the number of degrees of freedom of the particles. Here, we had d = 3. 

The same analysis can be conducted for the relativistic case, yielding the relativistic 
version of Eq. nil), 



d^v /(O ,v) j{v)mv^. 


(BIO) 


In the ultra-relativistic limit, we can set mv'^ mc^ in th e integrand and obtain the 
ultra-relativistic equation of state ( Landau fc LifshitO 19801) . 



(Bll) 
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